Gauge-invariant construction of quantum cosmology 
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We present and analyze a gauge-invariant quantum theory of the Friedmann-Robertson- Walker 
universe with dust. We construct the reduced phase space spanned by gauge-invariant quantities 
by using the so-called relational formalism at the classical level. The reduced phase space thereby 
obtained can be quantized in the same manner as an ordinary mechanical system. We carry out the 
quantization and obtain the Schrodinger equation. This quantization procedure realizes a possible 
resolution to the problem of time and observables in canonical quantum gravity. We analyze the 
classical initial singularity of the theory by evolving a wave packet backward in time and evaluating 
the expectation value of the scale factor. It is shown that the initial singularity of the Universe is 
avoided by the quantum gravitational effects. 

PACS numbers: 04.60.-m, 04.60.Ds, 98.80.Qc, 04.20. Fy 



I. INTRODUCTION 



On large scales, observational data show that the Uni- 
verse is well described by general relativity (GR). Since 
the Universe is expanding at present, it becomes smaller 
as it evolve backward in time and at last we encounter 
the initial singularity. However, it is thought that quan- 
tum gravitational effects become dominant at the Planck 
scale, where the Universe itself must be treated as a quan- 
tum object. One thus needs a quantum theory of the 
Universe in order to understand what happens near the 
initial singularity. Symmetry reduced models, e.g., the 
isotropic and homogeneous model P, 0, H, 0] ) have been 
canonically quantized by the Dirac quantization method, 
which leads to the Wheeler-DeWitt equation [1, 0] . Un- 
fortunately, this method causes the "problem of time" 
(see, e.g., [3] and references therein), that is, time evo- 
lution is lost in the following two senses. The first is 
the dynamics of the wave function The canonical 
formulation of GR leads to a constrained system and a 
Hamiltonian is of the form H — 'SiN^Ci, where Ci — 
are first-class constraints and iV* are arbitrary functions. 
Since the classical constraints are translated into the 
quantum constraints Ci^* = by the Dirac method, 
the Schrodinger equation becomes trivial for the phys- 
ical state: = = 0. Therefore, the dynamics 
of the wave function is lost in the quantum framework. 
The second is concerned with observables. In ordinary 
gauge theories, observables are gauge-invariant quanti- 
ties. Since a gauge- invariant quantity O is defined as hav- 
ing vanishing Poisson brackets with all constraint func- 
tions {Ci,0} — 0, the quantity O becomes a constant 
of motion by the Hamilton equation O = {H, O} — 0. 
Thus, the dynamics of observables is lost in both clas- 
sical and quantum theories if one restricts observables 
to gauge-invariant quantities. In fact, the problem of 
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what should be interpreted as observables in classical GR 
has long been discussed by many researchers, see for in- 
stance [I s [13 [n [11 [H [11. 

A possible way to overcome these problems is that one 
finds gauge-invariant quantities and a method to extract 
their physical evolution at the classical level, and then 
constructs a quantum theory based on them. The idea for 
constructing gau ge-invariant quantities which has been 
stressed in [9|, Il3| is that the relation between dynamical 
variables is gauge invariant even if they are gauge vari- 
ant, respectively. As the realization of the idea, a formal 
expression of gauge-invariant quantities in constrained 
systems was recently presented [isl . [l6j . This method is 
often called the relational formalism. Although the ap- 
plication of the formalism to full GR is hard in general, 
the situation is dramatically simplified if the constraints 
can be written in the so-called deparametrized form (see, 
e.g., [13] )■ If one applies the relational formalism to a 
deparametrized theory [l^, [l^, one can construct the 
reduced phase space coordinatized by gauge-invariant 
quantities and obtain a physical Hamiltonian which gen- 
erates the time evolution thereof. Then, one can quantize 
the reduced system in the same manner as in elementary 
quantum mechanics and obtain the Schrodinger equation 
because there are no constraints in the reduced phase 
space. 

In the present work, we construct and analyze a gauge- 
invariant quantum theory of the Friedmann-Robertson- 
Walker (FRW) universe without the problem of time. We 
consider the case when the matter involves dynamical 
dust coupled to gravity, where the dust action is intro- 
duced by Brown and Kuchaf [23| . The advantage of the 
dust is that one can deparametrize the Hamiltonian con- 
straint and extract a natural time variable which corre- 
sponds to the cosmological proper time when one solves 
equations of motion. Therefore, we can apply the rela- 
tional formalism to construct the reduced phase space 
of the FRW universe with dust, and then quantize the 
reduced system without dealing with the constraint. In 
order to investigate what happens near the classical ini- 
tial singularity, we find and analyze the solutions of the 
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constructed quantum theory of the Universe. We first 
construct a wave packet which allows semiclassical inter- 
pretation. Then, we evolve it backward in time by the 
Schrodinger equation and evaluate the expectation value 
of the scale factor. As a consequence, it is shown that 
the initial singularity is replaced by a big bounce by the 
quantum gravitational effectsD 

The organization of the paper is as follows. In Sec. 
im we review the relational formalism for deparametrized 
theories and present the classical theory of the FRW uni- 
verse with the dust in Sec. IIIII Then, we quantize the 
obtained system in Sec. llVi where we divide the section 
into three parts by the operator ordering of the Hamil- 
tonian. We analyze the dynamics of the Universe in Sec. 
IVl Sec. IVII is for the conclusion. 

In this paper, we adopt the following unit for the speed 
of light: c = 1. 



II. RELATIONAL FORMALISM FOR 
DEPARAMETRIZED THEORIES 

In this section, we summarize the relational formalism 
in the case where the constraint can be written in the 
deparametrized form. See [l^, [lB| for general cases and 
details. From now on, we assume that the system has 
only one constraint for simplicity. 

The key observation of the relational formalism to de- 
fine gauge- invariant quantities is as follows. Take two 
functions F and T on the phase space. Then, the value 
of _F at r = r is gauge invariant even if F and T them- 
selves are gauge variant. That is, one can interpret one 
of the functions T as a clock and consider the relation 
between T and other variables as time evolution. If we 
denote a phase space point hy x — {q°',Pa), the mathe- 
matical definition of the gauge- invariant quantity O'^p {x) 
as a phase space function is given by 



gauge orbit 



O'p^x) -.^ a%{F){x)\^t^(T)(x)-- 



(1) 



Here, denotes the action of the gauge transforma- 
tion generated by C, and t is a parameter along the 
gauge orbits. Thus, oi*(j{x) is the gauge flow generated 
by C starting from x. The action of on a func- 
tion is given by a\j{F){x) = F{a]j{x)), which is writ- 
ten as a series a*c(F)(a;) = Yln=Q ^{'^^P}{n){x), where 
{C,F}(o) := F and {C,f;}(„+i) := {C, {C, f^}(„)}. The 
definition ([1]) gives a manifestly gauge-invariant quantity 
because Op{x) is constant on each gauge orbit. Indeed, if 
X and x' are on the same gauge orbit, there is t' such that 



F 



0},{x') (see FigH]). 



T(a*^(i.'))=T 

A constraint equation C = is said to be of de- 
parametrized form if it is written as 




F=const. 



FIG. 1: The action of the map a\; on phase space points and 
the gauge- invariance of 0^p{x). 



with some phase space coordinates {q°' ,T;pa, P}. 
In the deparametrized theories, the reduced phase 
space is spanned by the gauge-invariant quantities 
{Oqa {x),Op^ (x)) associated with g° and Pa with the sim- 
ple symplectic structure 



{o^^4x),o;^{x)}^s^ 



(3) 



The physical Hamiltonian H is obtained by replacing 
and Pa in h{q°-,pa) with OJa(x) and Op^{x), 

H {Ol. ix),Ol (x)) h {0\. [x\ O;^ {x)) , (4) 

which generates the time evolution of the gauge- invariant 
quantities: 



dOUx) 
dr 



= {H,OUx)}. 



(5) 



As we have seen, one can construct the gauge-invariant 
quantities and extract their physical evolution in the de- 
parametrized theories by using the relational formalism. 



III. FRIEDMANN-ROBERTSON-WALKER 
UNIVERSE WITH DUST 

In this section, we shall formulate the FRW universe 
with dust. As was mentioned above, the advantage of 
the dust introduced by Brown and Kuchaf ^] is that 
one can deparametrize the system. 

The action of the dust is given by 



'5'dust — 



d^xV^pig'^'U^.U, + 1), (6) 



M 



C ^ P + h{q\pa) = 



(2) 



where the one-form is defined by = —{dT)^ + 
Wj{dZ^)j^ (j = 1,2,3), p is the rest mass density of the 
dust and g is the determinant of the metric tensor g^^, on 
the spacetime manifold M. The action is written by the 
variables 5^1,, p, T, and Wj. When the equations of 
motion hold, labels the flow lines of the dust particle 
and T is the proper time along them. By using the Dirac 
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algorithm for constrained systems, one can obtain the 
canonical form of the dust action 



consider the case of the flat universe, k 
is written in Cartesian coordinates as 



0. The metric 



'S'dust — 



), (7) 



with the dust degrees of freedom being T and and 
their conjugate momenta P and Pj, respectively, where 



Cdust = ^P^ + q-b J^^nstJ^dust^ 
D^^-^ ^PT^ + PjZ%. 



(8) 
(9) 



Here, N and N°-, respectively, are the lapse function and 
the shift vector and qab is the induced metric on a three- 
space E. It is worth noting that P is written by using the 
original variables as P = — p^det {qabjU^n^, where n^^ 
is the future-pointing unit normal to E. It is also shown 
that is a future-pointing unit timelike vector on shell, 
so that U^n^ < and the signs of P and p are the same. 
Recently, several applications of the relational formal- 
ism to systems with the dust have appeared [l9l . [2ll . [23 | , 
where the dust has negative energy in order for the physi- 
cal Hamiltonian to become positive definite. In our work, 
we choose the signs of P and p as positive in order to in- 
terpret the dust as a standard matter. Although this 
choice leads to a negative Hamiltonian, it is not surpris- 
ing because the contribution of the gravitational fields to 
the Hamiltonian is originally negative. 

In the case of the FRW universe, the dust action is cast 
into the form 



'dust 



dt / d^^T FT-iVCdust 



where 



a 



dust 



= p. 



(10) 



(11) 



Here, T corresponds to the cosmological proper time 
when one solves the equations of motion. Therefore, it 
can be interpreted as a natural time. This will be used 
as a clock in the discussion below. 

On the other hand, the action of gravity in the FRW 
spacetime is given by 



where 



'S'crra 



Cgra 



dt / d^X (pad ~ TVCgiav) : 



■A 

'I2a 



Aa3 

K 



3ka R 
K a 



(12) 



(13) 



Here, a is the scale factor, pa is its conjugate momentum 
defined as :— — A is the cosmological constant, 
R = Pradfl^ is a constant associated with the kinemati- 
cally incorporated radiation with energy density Piad a-nd 
the constant k = ±1,0 determines the curvature of the 
three-dimensional space. For simplicity, we shall only 



ds^ = -dt^ + a^{t)(dx^ + dy^ + dz^), 



(14) 



where we assume that the scale factor has a dimension 
of length and the coordinates are dimensionless. 

The total action of the FRW universe with the dust 
takes the form 

'S'tot — '5'grav + 'S'dust 

= J dtj d^x (pad + Pf- NCtot) , (15) 

C'tot — C'grav + Cdust — Cgrav + P- (16) 



Here, the constraint is of the almost deparametrized form 
except for the existence of the three-space integral. In 
general, since the volume of the three-space will diverge 
in the fiat case, we must somehow get rid of this diver- 
gence in order to deparametrize the system and consider 
the canonical quantization. 

We shall avoid the divergence and deparametrize the 
system by considering a compact universe. For simplic- 
ity, we only consider the case of three-dimensional torus, 
where we take a cube of coordinate range < a:, y, z < L, 
and identify the opposite faces. Then, the scale factor 
multiplied by L is the length of an edge of the cube, 
which, in fact, represents the physical size of the Uni- 
verse. If we denote the coordinate volume a.s V := 
J-^ d^x = L^, the change of the variables a' := V^a — aL, 

p'^ :— V^pa and P' :— VP absorbs the volume factor in 
the total action: 



Pa a 



CU^P' + hia',p',), 



P'f - NCL 



„/2 

'12a' 



Aa'3 R 



(17) 
(18) 
(19) 



Here, a' is the physical length of an edge of the Uni- 
verse (three-torus) and thus V^hys ■— o,'^ = a?L^ is the 
physical volume of the universe. So far, we have finished 
deparametrizing the system of the FRW universe with 
the dust. 

Let us apply the relational formalism to the system and 
obtain the reduced phase space coordinatized by gauge- 
invariant quantities. In the present case, T becomes the 
clock, that is, the evolution of all gauge- invariant quan- 
tities are measured by a relative relation with respect to 
T. As explained in SecHIl the reduced phase space is 
coordinatized by 



A{t) :=0:,(a;), Pa{t) :=OL(a;). 



(20) 



Recall that, the meaning of 01^, {x) is the value of a' at 
T = T. The symplectic structure is written as 



{A{t),Pa{t)} = 1. 



(21) 
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As can be seen from ([4]), the Hamiltonian which generates 
the time evolution of the gauge-invariant quantities is 
obtained by replacing a! and p'^ in h{a' ,p'g) with A{t) 
and Pa{t)- 



H = h{A, Pa) 



p2 

'12 A 



AA3 R 

1^ -J- 

K A 



(22) 



We have obtained the reduced phase space of the FRW 
universe with the dust. The system now is character- 
ized by (pij) and with no constraints, so that we 
can quantize the system in the same procedure as in the 
ordinary canonical quantization. Hereafter, we call the 
gauge-invariant quantity A{t) the scale factor. 



IV. QUANTIZATION 

Let us now quantize the system obtained in the previ- 
ous section. Our classical system is one-dimensional with 
the symplectic structure (1^ . In the ordinary procedure 
of the canonical quantization, the Poisson bracket is re- 
placed by the commutation relation among the operators 
corresponding to canonical variables: 



[A,pA]=ih. 



(23) 



We take the ordinary Schrodinger representation of the 
canonical commutation relation, that is, A acts on the 
state VE'(A) in the Hilbert space defined below by multi- 
plication, and Pa acts on 4'(A) by differentiation: 



i*(A) = A-^iA), Pa'^IA) 



(24) 



In general, there are many ambiguities in the choices of 
canonical variables, the Hilbert space and the operator 
ordering. As for the first one, we have chosen the most 
natural phase space variables, the scale factor and its 
conjugate momentum. Although this choice leads to neg- 
ative values of the scale factor by von Neumann's theo- 
rem, one can consider the restriction of the range of the 
scale factor A> in the quantum system. Under the re- 
striction. Pa fails to be self-adjoint in general. However, 
one can ensure the self-adjointness of the Hamiltonian by 
imposing boundary conditions on wave functions. Thus, 
we first define the operator ordering of the Hamiltonian 
(I22p and the measure on the space of wave functions, and 
then derive the boundary conditions on wave functions. 

We consider the three cases where is promoted to (a) 
Pa^Pa, (b) i (Pli + Pa^Pa + iPj) and (c) 



1 p2 
A^A- 



A. Operator ordering (a) — » Pa-jPa 

In the case (a), the Hamiltonian operator is given by 

K ^ 1 ^ AA^ R , , 

12 A K A 



and the Schrodinger equation = H'^ takes the form 
1 1 ai-X f KA^ R" 



ih—— = 
dr 12 



A dA^ A2 dA 



V 



A 



(26) 



We define the Hilbert space as TC 



dA) where 



IR+ represents the set of non-negative real numbers. With 
this inner product. Pa is Hermitian and in fact symmet- 
ric, but not self-adjoint. The condition for the Hamilto- 
nian to become Hermitian {^'i\H\'^2) (*i|-H^^|^'2> is 
satisfied when 



dA 



1 d^^y 



1 d^>-. 



I 



A dA^ dA 
1 d^^l 



dA 



A dA^ 



I 

A^~dA 



(27) 



Integrating the first terms on the both sides by parts, we 
obtain the equation 



1 d-^o 



A dA 



Ha A 



(28) 



If we assume the wave functions vanish at infinity, Eq. 
m yields ^a^L^o = 1§i*2L^o and the relation 
holds if each of the wave functions and ^2 satisfies 
the condition 



1 f rf^ 
VA [ ^dA 



= 0, 



(29) 



^=0 



where 7 is a real number. The Hamiltonian is indeed 
self-adjoint if its domain is restricted to wave functions 
which satisfy the boundary condition 

As we have seen, the Schrodinger equation and the 
boundary condition to ensure the self-adjointness of the 
Hamiltonian in the Hilbert space (R+ , dA) are derived. 
In Sec.V, we shall use the simplest boundary condition 
corresponding to 7 = 0: 



Va 



= 0. 



(30) 



A=0 



For later use, we present the dimensionless form of the 
Schrodinger equation. Since the characteristic length and 

time scales here are the Planck length / p = ^ ^ and the 

Planck time tp — \J we have the following dimension- 
less equation by using them: 



.d-^ 1 d^^ 
drj X dx^ 



1 / 
a;2 dx V 



Xx-^ 



xJ 



(31) 



where x 
V 27r / h- 



\2n) Ip 



hG 



A and r 
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B. Operator ordering (b) 



LA. 

A 



In the case (b), the Hamiltonian operator is given by 



H 



36 



Pi 



1 



1 



'A 



A 



A 



Ai3 R 



(32) 



and the Schrodinger equation is written in the form 
a* nh^ / 1 I 2 

ifl— = —— -r^—TTT TTTTr-r + TT-TT* 



dr 12 \AdA^ dA iA^ 
/A A3 R 



(33) 



Let us define the Hilbert space asH. — (M+ , dA) . Then 
Pa is Hermitian as well as in the case of ordering (a). 
The condition for the Hamiltonian ([5^ to be self-adjoint 
is also given by (P^ . 

The dimensionless Schrodinger equation is given by 

1 19* f 2 ^3 r\ . ^, 



a; Qa;^ dx 



ix^ 



C. Operator ordering (c) 



1 p2 



In the case (c), the Hamiltonian operator is given by 
= ^^P" + — + ^, (35) 



K 1 KA^ R 
12 A K A 



and the Schrodinger equation is written in the simplest 
form 

*V-T2Aa4^+(— 

Although we have defined the Hilbert space naturally as 
Ti. = L2(M_|_, dA) in the previous two cases, here one can- 
not obtain specific boundary conditions with that defini- 
tion. In order to avoid the difficulty, we shall choose the 
Hilbert space as H = , AdA) consisting of square 

integrable functions of A with respect to the measure 
AdA. With this inner product, Pa cannot be Hermitian. 
It can be shown that the boundary condition on the wave 
functions for H to be self-adjoint is ^|^^q = T^^I^^q' 
which is discussed in the context of the Wheeler-DeWitt 
theory [H^ [23] . The two simplest cases with 7 = and 
7 = 00 are 



(l)*L.o-0 (2) — 



0. 



(37) 



A=0 



The dimensionless Schrodinger equation is written as 
.d^ _ 

X dx^ 

will be dis- 



dr] 



1 / r\ ^ 

+ i^Xx^ -h - ) *. 



(38) 



The boudary conditions ([57)1 and Eq 
cussed in Appendix A. 



V. DYNAMICS OF THE UNIVERSE 

In this section, we shall analyze the dynamics of the 
Universe. We first consider the case (a) discussed in the 
previous section. We perform numerical calculations in 
order to solve the Schrodinger equation (pij) because one 
cannot analytically solve the equation. The procedure 
of the numerical calculations is as follows. We first set 
the initial wave function sharply peaked at some value 
of A. We next evolve it backward in time and evaluate 
the expectation value of the scale factor. The boundary 
condition considered here is pop . The numerical methods 
used here are the fourth-order Runge-Kutta method in 
the time integration and the midpoint difference method 
for the spatial differentiation. Although the equation to 
solve is (PT|) . it is not simple to control the errors in the 
whole range of x because of the rapidly growing potential 
term proportional to x^. Therefore, we change the spatial 
variable to y :— x^ , so that the Schrodinger equation ([31]) 
becomes 



drj 



9 52*(y,,7) 
4 



3 1 g^(2/,r?) 
Ay dy 



(39) 



We choose the initial wave function as 



^'(y,0)cxyexp(^-^^^-|^-zfco?/) , (40) 

which satisfies the boundary condition (l30|) . Although 
the other choices of the initial wave function are possible, 
one can obtain qualitatively the same results irrespective 
of them. Fig. 2(a) shows the comparison between the 



expectation values of the scale factor for the four cases 
(r,A) = (0,0), (0,1), (1,0), and (1,1), where we set the 
initial conditions as yo — 7, ko = 20, and a = 0.5. The 
absolute values of the wave function when r = A = 1 
are plotted in Fig. 



3(a) 



The effect from the radiation 
is so small that the plots for the models with or without 
the radiation almost completely overlap. Therefore, we 
only show the two cases when the matter involves the 
radiation, with or without the cosmological constant, in 
Fig. 



2(a) 



As the figure indicates, the expectation value 
of the scale factor has a nonzero minimum, that is, the 
initial singularity is replaced by a big bounce. It also can 
be seen from Fig. 2(a) that the effect of a cosmological 



constant slightly decelerates the expansion of the Uni- 
verse. We note that the values of the cosmological con- 
stant and the energy of the radiation used here are much 
larger than those suggested from the astrophysical ob- 
servations. Therefore, the actual contributions from the 
cosmological constant and the radiation are still smaller, 
so that the results obtained here can be thought as the 
universal feature. 

We also perform the above calculations in cases (b) 
and (c). The expectation value of the scale factor and 
the absolute value of the wave function in case (b) are 
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FIG. 2: The expectation value of the scale factor {A) as a 
function of time r for the four cases (r, A) = (0, 0), (0, 1), (1, 0), 
and (1,1), where r represents the dimensionless parameter 
associated to the total energy of the radiation while A is the 
dimensionless cosmological constant. We show only the cases 
when the matter involves the radiation because the plots in 
the models with or without the radiation almost completely 
overlap. The solid line represents the cases with the cosmo- 
logical constant, while the dashed line is for the cases without 
the cosmological constant. Figs. (a), (b), and (c) correspond 
to the three kinds of the operator ordering discussed in Sec lIVI 




(a) 





(c) 



FIG. 3: The absolute value of the wave function is plotted as 
a function of the time r and the scale factor A in the case 
when the matter consists of the dust, the radiation, and the 
cosmological constant. Figs. (a), (b), and (c) correspond to 
the three kinds of the operator ordering discussed in Sec lIVI 



illust rate d in F igs. 2(b) and 3(b) respectively, and Figs. 
2(c) and 3(c) are for case {cf. As Fig. [5] indicates, we 
obtain almost identical behavior, qualitatively and quan- 
titatively, of the expectation value of the scale factor ir- 
respective of the operator ordering. We see from Fig. [3] 
that the peak of the wave packet in case (c) is lower in 
comparison than with cases (a) and (b). 



A special feature of case (c) is that one can analytically 
solve the Schrodinger equation, construct wave packets 
and compute the expectation value of the scale factor 
when the matter consists only of the dust. We show 
that in Appendix A and the same calculations appear 
in the context of the Wheeler-DeWitt theory [23| , where 
observables are not gauge-invariant. 



VI. CONCLUSION 

We have constructed and analyzed a gauge-invariant 
quantum theory of the FRW universe with dust. In or- 
der to obtain the quantum theory free from the problem 
of time and observables, we have used the following two 
properties: (i) the Lagrangian with the dust introduced 
by Brown and Kuchaf enables one to deparametrize the 
system, that is, the Hamiltonian constraint has the form 
C = P + h{q'^,pa) = and (ii) deparametrized theo- 
ries allow the construction of the reduced phase space 
coordinatized by gauge-invariant quantities through the 
relational formalism. We have made a natural choice 
of the time variable which agrees with the cosmologi- 
cal proper time when one solves the equations of mo- 
tion. Then, we have quantized the reduced system and 
obtained the Schrodinger equation. In order to analyze 
whether the initial singularity is avoided by the quantum 
gravitational effects, we have constructed wave packets, 
evolved them backward in time and evaluated the expec- 
tation value of the scale factor. It has been shown that 
the expectation value of the scale factor never goes to 
zero, that is, the initial singularity of the Universe is re- 
placed by a big bounce. We have obtained qualitatively 
and quantitatively almost identical results irrespective of 
three operator orderings. Although we have used a spe- 
cific boundary condition in Sec.V, it is not necessarily re- 
sponsible for the results. Indeed, it is shown in Appendix 
A that a bouncing universe is obtained under two differ- 
ent boundary conditions. The construction of a quantum 
theory discussed in this paper can be extended to other 
models, e.g., the FRW universe with a scalar field, and 
we shall study these problems in future works. 



APPENDIX A: DUST UNIVERSE IN THE CASE 

(C) 

In case (c), when the matter consists only of the dust, 
one can analytically solve the Schrodinger equation, con- 
struct wave packets and compute the expectation value 
of the scale factor. The boundary conditions considered 
here are (l37l) . From the Schrodinger equation pS)) . the 
energy eigenvalue equation reads 



1 d'^ipix) 



X dx^ 

The general solution of the equation (lAip is 
V'(z) = CiAi(-z)-hC2Bi(-0), 



(Al) 



(A2) 
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where z := {—E)^x, Ai and Bi are the Airy functions and 
Cj's are arbitrary constants. If z has the sign opposite 
to a;, that is, if £'> 0, it is easy to see from the behavior 
of the Airy functions that there are no wave functions 
which has a finite norm and satisfy the boundary condi- 
tions ([37|) . Therefore, the energy eigenvalue E must be 
restricted to the range E < ^. This fact is reasonable 
because the total energy, which consists of the gravita- 
tional energy and the matter energy, is constrained to 
vanish and E can be interpreted as the energy of grav- 
ity. Under the restriction < 0, we can transform the 
wave function (|A2[) by using the relation between the 
Airy functions and the Bessel functions J^\: 



^{z) = 



(A3) 



It follows from the behavior of the Bcssel functions for 
small argument that the allowed wave functions under 
each of the boundary conditions (|37|) are 



(A4) 



(1) ^«(z) = ViJi ( 

(2) i,%\z) = ^zJ_.^{^ 



2 3 

3 



Although these solutions are not normalizable with re- 
spect to the measure AdA, we can construct wave packets 
which have a finite norm by superposing them. This sit- 
uation is similar to the case of a free particle in quantum 
mechanics where the energy eigenstates are not normal- 
izable and one often constructs a wave packet to see the 
motion of the particle. 

The general solutions of the Schrodinger equation 
jO* _ .^ypj^tgQ [t^ tj^g form 

on X ax'^ 



where 



e = ^^/—E and I = 1,2. If we choose 

the functions C'^^He) to be C'^^\e) = e""'' and 

C"*-^-' = £"3 6""*^ where a is an arbitrary positive 

constant, we can use the formula [25] for the Bessel 

function e^"-'^^ x"^^ J^{bx)dx = f^^^y+i e^^ ,Re{a) > 
0,Re(i^) > — 1. Then, the wave functions take the form 
of the following wave packets: 



¥^\x,ri) 



X _x3 

e ^, 



-e , 



(A6) 



where /3 := a — i^rj. We interpret the wave packets 
as the states of the Universe. 

The expectation value of x is given by (z)^^) = 
Iq° I'^^^^l'^x'^dx/ l^f^^^pxiix and we can calculate this 
integration with respect to the wave packets: 



(x) 



(1) 



2r 



3r 



2a 



3.2§r(f) 



2a 



(A7) 



These solutions show that the expectation value of the 
scale factor never goes to zero, that is, the initial singular- 
ity of the Universe is avoided by quantum gravitational 
effects. The asymptotic behavior of the expectation val- 
ues for large x becomes 



4'(^)(x,77) 



C^-^1{E)e~'^'^'ij4^dE 



C"(i)(g)e»f^^£|+ij^, (ext) de, (A5) 



(A8) 



Thus, they are in good agreement with the classical tra- 
jectories when the Universe is sufficiently large. 
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